Squeezing of electromagnetic field in a cavity by electrons in Trojan states 
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The notion of the Trojan state of a Rydberg electron, introduced by I.Bialynicki-Birula, 
M.Kaliriski, and J.H.Eberly (Phys. Rev. Lett. 73, 1777 (1994)) is extended to the case of the 
electromagnetic field quantized in a cavity. The shape of the electronic wave packet describing the 
Trojan state is practically the same as in the previously studied externally driven system. The fluc- 
tuations of the quantized electromagnetic field around its classical value exhibit strong squeezing. 
The emergence of Trojan states in the cylindrically symmetrical system is attributed to spontaneous 
symmetry braking. 



I. INTRODUCTION 

The possibility of creating stationary, nondispersive, 
localized, wave packets describing a highly excited elec- 
tron in a hydrogen atom strongly driven by circularly po- 
larized microwave radiation was predicted theoretically 
severalyears ago Q and confirmed in numerous publica- 
tions |g-0] (for recent reviews of the subject see 
Such electronic states are called Trojan wave packets by 
analogy with the cloud of Trojan asteroid in the Sun- 
Jupiter system. 

In all previous studies the microwave field was treated 
as an external, classical wave. Dressing of an electron by 
such a wave of a suitably chosen intensity and the fre- 
quency equal to the Kepler frequency of the electron on 
the Rydberg orbit makes the Trojan wave packets highly 
stable. Their life-time is of the order of one second P,|i"o|, 
which makes them an interesting object of study for the- 
oretical and perhaps even for practical reasons. 

In the present paper a similar problem of nondisper- 
sive electronic wave packets is studied for an atom in- 
teracting with the quantized electromagnetic field. Such 
an approach allows for a fully dynamical treatment of an 
autonomous atom-field system. It automatically includes 
a back reaction of the atom on the electromagnetic field. 
Thus, one can study both the dynamics and the statis- 
tical properties of the electromagnetic radiation. Our 
study fully confirms the existence of Trojan states of the 
Rydberg electron in this new regime with almost exactly 
the same shape of the wave packet. The back reaction 
of the electron on the electromagnetic field pushes the 
field frequency off resonance. The quantum fluctuations 
of the electromagnetic field exhibit strong squeezing. 



II. HYDROGEN ATOM IN A CAVITY 

Anticipating the role of highly populated, discrete 
modes of the microwave field in the formation of Tro- 
jan electronic states, we consider a hydrogen atom in a 
microwave cavity. In the presence of a cavity we can 
separate a finite number of relevant degrees of freedom 
whereas in free space we would have to deal with a contin- 
uous spectrum which precludes the existence of localized 
stationary states of the system. 

To allow for the rotational symmetry of the atom-field 
system we choose a cylindrical cavity. Its dimensions will 
be large enough to justify the dipole approximation in the 
coupling of hydrogen atom with the lowest cavity modes. 
The atom placed in the middle of cavity interacts only 
with TEi„ modes. For definitness we choose the two (de- 
generate) lowest modes of this type (n = 1) (labeled by 
X and Y) for which the mode functions have the form 

E x = iA/w sin ^— e z x Vj_Ji (xur/R) sin <p, (la) 



B x = — cos— V_LJi(a;iir/i?)sin^, (lb) 

Ii T 

E r = — iNuj sin — e 2 x V^Ji (xur/R) cos (p, (lc) 



B = — — cos — Vj_Ji(xnr/i?) cos^, (Id) 
T T 

where R and L are the radius and the length of the cavity 
and x\i is the first (the smallest) solution of the equa- 
tion dJi(x)/dx — 0. The z-axis is taken along the cylin- 
der axis and Vj_ = (d/dx,d/dy). The frequency of the 
modes and the effective wave vector are given as 
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u = j i (x 2 11 + (nR/L) 2 fz, 
k = ^{lu/cY + (n/L)2. 



(2) 
(3) 



The value of the normalization constant J\f has been ob- 
tained in Ref. pi, 



TV = 



k 2 R 2 V 2tt€ Llu(1 - l/x? x )J?(a:ii) ' 



(4) 



from the requirement that the energy per one photon in 
a mode is equal to Tiuj. 

At the position r.4 = (0, 0, L/2) of the center of the 
atom the orthogonal field vectors E x and E F point re- 
spectively in x and y direction and are given by simple 
formulas: 



E^(r A ) = 
E r (r A ) = 



-iMwxw 

2R 
-iJ\fu>xn 



2R y ' 
B x (r A ) = 0, B Y (r A )=Q. 



(5a) 

(5b) 
(5c) 



The relevant part of the electric and magnetic field in the 
cavity can be written in the form 

E = -E x a x + E x *a* x + E Y a Y + ~E Y *a* Y , (6a) 
B = B x a x + B x *a* x +B Y a Y + B Y *a Y , (6b) 

where a x and a Y are the dimensionless mode expansion 
amplitudes. 

In the laboratory frame the dynamics of the atom-field 
system is governed by the Hamiltonian 

2 2 
H L = ^---^ er-E(r A ) 

+ 1 J (e E 2 + B 2 / Mo )d 3 r, (7) 

where r = [x, y, z) is the position of the electron rela- 
tive to the center of the atom r^. The Hamiltonian Hl 
describes the mutual interaction of the atomic electron 
with the chosen cavity modes. We can rewrite Hl using 
the amplitudes a x and a Y , or more conveniently, using 
their real combinations: 

— z — 2 

Px = ~^(a x -a* x ), P y = -j=(a Y ~ a Y ), (8a) 



Q x = -j={a x + a* x ), Q y = -y=(a Y + a Y ), (8b) 

where the dimensionless vectors P and Q represent the 
electric field and the magnetic induction, 



Hr = 



p 2 e 2 



e £r.P + ^(P 2 + Q 2 ). (9) 

2m 47reo^ 2 



The field amplitude £ is 



£ = 



Afuix 
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RV2 



(10) 



We have found it convenient to use the natural units for 
our problem derived from the field f requen cy for the en- 
ergy, length, and momentum: Tiuj, ^Jh/muj, \fTimuj. The 
Hamiltonian (^) in these units takes on the following form 

2 + Q 2 



ifc-*--J- 7 ,r + i a 

where the dimensionless parameters q and 7 characteriz- 
ing the strength of the Coulomb field and the atom-field 
coupling are 



rtiLO 

~h~' 



e£_ 

Tiuj 



(11) 



III. CLASSICAL SOLUTIONS 

We are interested in special solutions corresponding to 
the Trojan states in the external electromagnetic wave 
introduced in Ref. (lj. Since these states describe elec- 
tronic wave packets rotating around the nucleus along 
circular orbits, we transform the Hamiltonian Hl to the 
frame rotating around the z-axis with the angular veloc- 
ity f2. The transformed Hamiltonian is 



H 



2 



71- P 



Q 2 



k(M* + M*), (12) 



where k = Cl/ui. The z-components of the angular mo- 
menta of the electron and of the electromagnetic field 
are M Z A = xp y - yp x and Mj = (Q x P y - Q V P X )- In 
this frame, the rotational states will appear as station- 
ary states of the Hamiltonian. We would like to stress 
that the Hamiltonian H cannot be identified with the 
energy because of the appearance of the inertial forces in 
the rotating frame. 

To emphasize the rotational symmetry of our problem 
we introduce following variables for the electromagnetic 
field 



Pa- 



Qx Py 

V2 ' 

Qy H~ Px 

~V2~' 



P_ 



Qx H~ Py 

~V2~ 



Px Qy 



V2 ' 



(13a) 



(13b) 



corresponding to the left and right circular polarization. 
In terms of these variables the Hamiltonian takes on the 
form 



p q 7 

2 r V2 

1 + K 



H=^--- -h=(x{P+ + P_) + »(Q_ - Q+)) 



(14) 
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The kinetic part of the field Hamiltonian is made up of 
two terms: co-rotating and counter-rotating. Linear sta- 
bility analysis shows that both parts are necessary for the 
existence of the nontrivial equilibrium solution. From the 
Hamiltonian (|12|) we derive the evolution equations: 



x = p x + ny, 



y =p y - KX, 



Z = Pz 



Q+ = (1 + k)P+ - fx/ V2, 
Q_ = (1 - «)P_ - fx/V2, 



qx j(P+ + P-) 
p x = -^3 + ^= + Kft,, 



. _ _ qy , i{Q- - Q+) _ 

Py — ~~ 3 1 7^ Kp x , 



V2 



qz 



P + = -(l + K )Q+- 7 y/V2, 



-(l-K)Q-+",y/y/2. 



(15a) 
(15b) 
(15c) 
(15d) 
(15e) 

(15f) 

(15g) 

(15h) 

(15i) 
(15j) 



The time-independent solutions of these equations de- 
scribe the stationary states of our system. Equating the 
left hand side of Eqs. (|l5| ) to zero, we obtain 



ro cos tp, y eq — r Q sintp, z eq = 0, 



(16a) 



p eq = — nr sirup, p e y q — nr cos tp, p eq — Q, (16b) 



_ -77-0 sin tp eq _ -77-0 sin ip 



(k + 1)V2' 



PI 



jro cos tp 
(k + l)%/2 



peg 



(k- 1)V2' 



— 7fo cos 93 
(k-1)V2' 



(16c) 



(16d) 



In addition, the equations (151) and (|l5gj) give the equi- 
librium condition: 



1 _ 2 

o — AC 

7'° 



7 



K 2 - 1 



(17) 



This equilibrium condition can be used to express the 
equilibrium radius ro in terms of the frequency of the 
cavity mode u> and the frequency of rotation O 



ro(fi) 



£1/3 



7 



-1/3 



(18) 



(fi/w) 2 - 1, 

or, alternatively, to express the frequency of rotation in 
terms of to and ro. The equilibrium condition ([17] 
two solutions for f2, denoted by fP(ro) and Jl < (ro) 



has 



n>(r ) = ^1 + q/rl + -/(W/r 3 )2 + 4 7 2 , (19a) 



^<(r ) = ^\f+ <l/r 3 o - ^0~q/r 3 o) 2 + 4 7 2 - (19b) 

Both solutions exist for all values of ro. The solu- 
tion fP(ro) (fP(ro)) gives the frequency that is always 
higher (lower) than the cavity frequency u (see Fig. [I]). 
The higher frequency (larger centrifugal force) requires 
the electric field to be directed towards the nucleus, 
whereas the lower frequency requires the field pointing 
outwards. The first case corresponds to the Trojan states 
whereas the second to the so called anti- Trojan states. In 
the previous study [|| , when the electromagnetic field has 
been treated as a given external wave, the anti- Trojan 
states were found to be classically unstable. The classi- 
cal stability obtained in the present study is due to the 
detuning from the exact resonance. Since the localization 
of the electronic wave packet is much worse for the anti- 
Trojan states (classically, the trajectories in the rotating 
frame are spread almost evenly around the whole circle, 
cf. Fig. ^j, we will restrict ourselves to the Trojan states 
only. Hence, in what follows, we shall only consider the 
solution fP(ro). 

Note, that for u> = (i.e. k = 1), we have only a trivial 
result r eq = p eq = P eq = Q eq = which in the classical 
model of an atom, means that "the electron has fallen 
onto the nucleus" and the electric field is zero. Thus, 
every nontrivial solution requires the presence of a de- 
tuning (lu Q) between the cavity frequency and the 
Kepler frequency. This phenomenon is known as the fre- 
quency pushing and is a direct consequence of the mutual 
atom-field interaction. This detuning has been absent in 
all previous approaches where the atom was driven by an 
external wave. 

The equations ( p"5| ) have a continuum of time- 
independent solutions that can be labeled by ro and the 
angle tp in the x-y plane. These solutions describe in the 
laboratory frame a classical electron circulating around 
the nucleus at the distance ro- The orbit of the electron 
is confined to the x-y plane. The electron is dressed by 
the classical electromagnetic field 



E 



£y/2 e£ro 



1 huj 



(sin tp e y + cos tp e x ), 



(20) 



which has a resonance dependence on the parameter k. 
Note that the electric field changes its sign when the fre- 
quency of rotation passes through the resonance. 

Next, we expand the Hamiltonian around a time- 
independent solution and investigate its linear stability. 
The motion will be stable if all eigenfrcquencies are real. 
The characteristic equation for this problem has the form 



3 



A 2 (A 2 -g r )(A 6 -(4 + , Zj , + 2)r)A 4 

+ (5 - 3/5g r + q 2 /2 + 4 7 2 + (4 + 5/2g r r)A 2 



-(2 + 5/2q r - 5q 2 +q 3 r + 8 7 2 + 14g r7 2 + (2 + 7/2q r + q 2 /2 + 8 7 2 )r))) 



= 0, 



(21) 



where q r = q/r^ and T = y/l — 2q + q 2 + 4 7 2 . The first 
(A = 0) frequency in our problem corresponds to the 
rotation of the whole system and it is a reflection of the 
rotational symmetry. The second frequency (A = \J q/r^) 
corresponds to the motion in the z-direction that (in the 
linear approximation) is decoupled from the motion in 
the x-y plane. The remaining 3 frequencies correspond to 
the motion of the electron coupled to the electromagnetic 
field. We shall not produce the analytical expressions for 
these eigenfrequencies but in Fig. [| we plot the region of 
stability in the R-tq plane. 

The stability can also be studied numerically and the 
calculations of the classical trajectories fully confirm 
the stability of the equilibrium solution. In Fig. ^ we 
plotted the projection of a typical electron trajectory 
on the x-y plane for the time interval (1400 T, 1500 T), 
where T = 1/uj. The trajectory started at the equilib- 
rium position r = ro (1,0,0) with the initial momenta 
p = mcjro(0.02, n + 0.07,0.02). As we see, the electron 
follows a rather complicated, but bounded, trajectory. 
Obviously, if we choose p(t — 0) sufficiently large the 
electron will eventually leave the vicinity of the equilib- 
rium point. In Fig. || we show the z-p z cross-section of 
the phase space for the same trajectory. This phase-space 
trajectory resembles the trajectory of a simple harmonic 
oscillator. Indeed, as we have seen, in the linearized evo- 
lution equations, the motion in the z-direction is purely 
harmonic. Thus, the interesting dynamics of the elec- 
tron is found in the motion confined to the x-y plane 
and in what follows we shall treat our problem as two- 
dimensional. 

Since our system is conservative, it has a well defined 
energy H^. We have calculated its value E(k) for all 
those solutions that in the rotating frame arc determined 
by Eqs. (|l6|). This energy is given by the formula 



mu ?q(k) / 7 



r(5K 2 -3) — k* 



(22) 



(k 2 - l) 2 

and is plotted in Fig. |as a function of k. The infinite 
growth of the energy near the resonance (k = 1) expresses 
the phenomenon of the frequency pushing. 



IV. QUANTUM EFFECTS 

In order to study the quantum effects for the electron 
as well as for the electromagnetic field we will apply the 
procedure of the quantization around the classical solu- 
tion (16). A similar quantization method has been used 
before, for example in nonlinear optics to describe quan- 
tum fluctuations around the classical solitons in fibers 



p4| . Here, the quantization will lead to the description 
of the electron in terms of a quantum mechanical wave 
packet orbiting along the classical trajectory and, at the 
same time, will reveal quantum fluctuations of the elec- 
tromagnetic field around its classical value. 

As a starting point we choose the Hamiltonian ( [Ti] ) in 
which all variables are treated as operators and we ex- 
press them as sums of their classical parts and the quan- 
tum corrections f = r eq + r, p = p eq + p, Q = Q 9 + Q, 

eq 

P = P + P. The classical parts represent equilibrium 
solutions ( |l6| ) found in the preceding Section. In order to 
simplify the notation, we have not attached any labels to 
the operators of quantum corrections (r, p, Q, P). Next, 
we expand the Hamiltonian around the classical equilib- 
rium solution neglecting all terms higher then quadratic 
in the quantum corrections. To proceed along these lines, 
we have to choose one solution, labeled by c/?o, from the 
whole family of equilibrium solutions. Making this choice 
we break the rotational symmetry of the Hamiltonian. 

This mechanism of selection of a specific classical 
solution resembles the spontaneous symmetry break- 
ing. Spontaneous symmetry breaking is present in many 
branches of physics. It explains the appearance of de- 
formed nuclei, the formation of magnetic domains in fer- 
romagnetic materials, or the emergence of Higgs particles 
in the Glashow-Weinberg-Salam model of electroweak in- 
teractions. In all these cases the symmetry is broken 
by the choice of a particular ground state. In our case, 
however, we do not break the symmetry by choosing a 
ground state but by choosing an equilibrium state of the 
Hamiltonian that is very far from the ground state of the 
system. 

Once we have chosen some , we can rotate the frame 
of reference, so that the direction given by ifQ is along the 
x axis. The quadratic Hamiltonian is 



H Q = \- fy x ( P + + P -) + - 3+) 



(P 2 + Q 2 _) 



2 9 2 2 

2 2 Q K V qK Z 

qn x -\ 1 K(xp y -xp x ), (23) 



where the parameter q, the ratio of the Coulomb force to 
the centrifugal force, 



AireomrQn 2 TqK 2 ' 



(24) 



has been introduced to achieve the full correspondence 
with the notation used before in the description 
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of Trojan states. Note, that in this Hamiltonian the 
quadratic term qn 2 x 2 enters with the negative coefficient. 
If it were not for the rotational term, such a Hamiltonian 
would not have any stable points. In our case, however, 
the stability can be achieved for a particular choice of 7, 
q and k . 

We look for the fundamental solution of the 
Schrodinger equation with the Hamiltonian Hq in the 
form of a four-dimensional Gaussian function 

i/> = Nexp(~'X.-A-x\, (25) 



— 2k 2 q — a 2 x + 2«ai2 + a 2 2 — 27013 + o 2 3 (l + k) — 27014 + a 2 4 (l — k) = 0, (27a) 

a n a 12 + a 12 a 2 2 ~ 7023 + oi 3 a 2 3 - 7 a 24 + ai4«24 + « (-an + o 22 + 013^23 - a 14 a 24 ) = 0, (27b) 

011013 + ai2»23 - 7033 + 0,130,33 - 7034 + ai4«34 + K («23 + ai3«33 - a 14 a 34 ) = 0, (27c) 

anou + ai 2 a 2 4 - 7034 + 013034 - 7044 + ai 4 a 44 + k (a 2 4 + ai3«34 - 014044) = 0, (27d) 

n 2 q — 2nai2 + a 2 2 — o 22 — 023(1 + «) — a 2 4 (l — k) = 0, (27e) 

- 7 - «ai 3 + ai 2 ai 3 - a 22 o 2 3 - 023033 - ^023033 - a 2 4a 34 + ^024034 = 0, (27f) 

7 - na 14 + ai 2 ai4 - a 2 20 2 4 - 023034 - ^023034 - a 2 40 44 + Ka 2 40 44 = 0, (27g) 

1 + k + a 2 3 - a 2 3 - (4(1 + k) - a 2 34 (l - k) = 0, (27h) 

ai30i4 - 023024 - a 34 ((1 + /«) a 33 - (-1 + k) a 44 ) = 0, (27i) 

1 - k + a\ 4 - a\ 4 - a 2 44 (l - k) - o 2 4 (l + k) = 0. (27j) 



We can easily solve these equations numerically, but 
first we want to find a perturbative solution. In order 
to do that we write the coupling c onstant in the form 
7 = 7\/k — 1. Obviously 7 = k^J (1 — q)(n + 1), and 
we will treat 7 as a small parameter. Typical values of 
the parameters are k = 1.0000001, q — 0.95625 which 
give 7 = 0.06. One can ask why we can not treat 7 (or 
even simpler, n — 1) as a perturbation parameter. How- 
ever, if we do so we face a problem: the coefficients of 
the perturbation series are growing, since they behave 
as 1/ V 'k — 1. When we tend with k — 1 to zero we hit 
exactly the resonance point and the perturbation expan- 
sion becomes meaningless. On the other hand, when 7 
is chosen as an expansion parameter, all large contribu- 
tions to the coefficients in the perturbation expansion 

a ij = a if + 7 a ij^ + l 2a if + • • • cancel out. 

We calculated analytically the coefficients up to the 



X = (x,y,Q+,Q-), and 



I an ia 12 ia 13 ia u ^ 

*Oi2 022 O23 024 

ia 13 a 23 a 33 a 34 

V »oi4 a 2 3 a 34 GS44 / 



(26) 



Inserting this ansatz into the Schrodinger equation wc 
obtain 10 algebraic, nonlinear equations for the parame- 
ters an: 



second order but we present here the analytic formulas 
only in the zeroth order and numerical values of the first 
and the second order corrections. 



(28a) 
(28b) 
(28c) 



(0) _ / (l + 2g)(4g-9g + 8-8s(g)) 
a n ^2 - 

(0) = 2 + q- 25(g) 

12 3g 

(o) _ / (l-g)(4g-9g + 8-8s(g)) 

°22 -«y ^ . 

where s(q) = ^1 + q — 2g, 

= 0, 4? = 0, ai? = 0, flM = 0, 
4°3 } - 1. «i°4 } - 0, = 1. (29) 
Thus, in the zeroth order, the electronic part of the wave 
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packet is exactly the same as in the case of externally 
driven Trojan wave packet El. The electromagnetic part 
has a form of a coherent (nonsqueezed) state. 

Higher corrections are due to the mutual interaction 
between the field and the atom. Numerical values of the 
parameters et^ are calculated for the cavity parameters 
L = 1 cm, B = 0.32 cm, which give u> — 197 GHz and 
7 = 3.24 x 10 -7 . The detuning k is chosen in such a way, 
that the value of q is optimal, q = 0.95625. As shown in 
Ref. §, the wave packet is then maximally concentrated 
around the equilibrium point and its center is located 
at ro = 3600ao (ao is the atom Bohr radius). The ex- 
pansion coefficients calculated up to the second order are 
presented in Table 1 . In this order we observe the effect of 
the back reaction of the electron on the electromagnetic 
field. However, the coefficients an,ai2, and 022 charac- 
terizing the shape of the electronic wave packet are the 
same as in the zeroth order within the assumed accuracy. 

Finally, we present in Table 2 the results of a direct 
numerical solution of our set of equations. As we see, 
almost all the coefficients have been obtained correctly 
already in the second order of perturbation theory. Only 
the 044 differs from the exact numerical solution. This 
can be attributed to the very slowly convergent perturba- 
tion series for this particular coefficient. The coefficients 
cti3, C23, ai4 and (Z24 describing the mixing of the atomic 
part of the wave packet with the field part are not zero. 
Moreover, the electromagnetic field is strongly affected 
by the interaction with the atom; the coefficient 044 is 
significantly different from its value in the zeroth order. 

The four-dimensional Gaussian wave packet (25) with 



the coefficients calculated numerically describes the 
fundamental state of the mutually interacting atom-field 
system. Owing to its Gaussian form, this state satu- 
rates the multidimensional generalized uncertainty rela- 
tions for the complete atom-field system (see Q). The 
smallness of the coefficients 013, 023, ai4 and 024 expresses 
the fact that the field and the atom are only very weakly 
correlated in this state. As a result of this, the saturation 
of the uncertainty relations is almost exact, separately for 
both parts of the wave function. The average values of 
second moments the electronic variables calculated with 
the numerical values of the coefficients taken from the 
Table 2 are given in Table 3. This Table exhibits the 
existence of correlations between the variables in the x 
and the y directions. This requires the use of generalized 
uncertainty relations for a two-dimensional system in the 
form (cf. Ref. |]|) 

(xx) (pxPx) + (xy) ipxPy) 

1 n 2 

~^{{xPx +Pxx) 2 + (xp y +p y x)(p x y + ypx)) > y, (30) 

(yy) (pyPy) + (xy) {pxP v ) 

1 h 2 
-^((ypy+Pyy) 2 + {ypx + p x y){p y x + xp y )) > —. (31) 

Upon substituting the values taken from the Table 3 we 
find an almost exact saturation of these relations. 



Now, we turn to the description of the quantum corre- 
lations for the electromagnetic field in our fundamental 
state of the atom-field system. The second moments for 
the field variables are given in the Table 4. These values 
of the correlations imply an almost complete decoupling 
between the co-rotating and counter-rotating modes so 
that the uncertainty is almost saturated separately for 
each mode. The state of the field in the counter-rotating 
mode is a coherent state but the state in the co-rotating 
mode is highly squeezed; the ratio of the correlation for 
the two quadratures Q_ and P_ is about 3.5 x 10 4 . How- 
ever, the fluctuations of the field are still small as com- 
pared to the field value Pt q = 1.5 x 10 6 . The plots 
of the Wigner function in Figs. [?] and || illustrates the 
difference in quantum fluctuations between the counter- 
rotating and the co-rotating modes. 



V. DISCUSSION 



We have shown that the dynamical treatment of the 
relevant modes of the electromagnetic field enables one 
to reproduce exactly the properties of the Trojan states 
studied previously in the presence of a given, external 
wave. However, there appear new features totally absent 
in the previous studies. First, the back reaction of the 
electron on the electromagnetic field causes the detun- 
ing from the exact resonance. As a result, the stability 
region covers now also the anti- Trojan states, that were 
before found to be classically unstable. Second, our anal- 
ysis has shown that to achieve the equilibrium state of 
the mutually interacting atom-field system we must take 
into account both polarization modes of the field: co- 
rotating and counter-rotating. The inclusion of only the 
co-rotating mode, as proposed in Ref. JT3], is not suffi- 
cient to achieve an equilibrium state. 

The quantization procedure adopted by us in this work 
consists in quantizing the corrections to the classical so- 
lution. These quantum corrections describe the shape of 
the electronic wave packet and the quantum fluctuations 
of the electromagnetic field around its classical value. 
The field fluctuations turn out to be significantly differ- 
ent for the two modes: for the counter-rotating mode the 
fluctuations are as for the vacuum state whereas for the 
co-rotating mode they exhibit strong squeezing. 

The choice of one particular solution from the class of 
equivalent classical solutions spontaneously breaks the 
rotational symmetry of the initial Hamiltonian. The 
method of quantization around the classical solution used 
here can be also applied to a similar problem of electronic 
Trojan states in a polar molecule ]l3[| . In this case, the 
role of the electromagnetic field is played by the rotating 
molecular dipole. The application of our method would 
require the dynamical treatment of the relevant molecu- 
lar degrees of freedom. 
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TABLE II. Coefficients characterizing the fundamental 
state of the atom-field system calculated numerically. 



an = 0.51160 0,12 

033 = 1 034 

ai 3 = 7.50 x 10~ 7 aw 

023 = —5.33 x 10~ 7 024 



0.78164 
1.40 x 10~ 12 
4.668 x 10~ 6 
-1.34 x 10~ 6 



a 22 = 0.06 2 70 
a 44 = 0.00532 
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TABLE III. Correlations of positions and momenta for the 
electronic variables calculated in the fundamental state of the 
atom-field system. The position variables are measured in 
units of the electron orbit radius ro and the momenta are 
measured in the corresponding unit mQro. 



(xx) = 0.01595 
(yy) = 0.13014 
(p x p x ) = 0.08369 
{PyPy) = 0.01026 

{xy) = 



{xp x + p x x) = 

(VPy + PyV) = 
(xpy+p y x) = -0.02493 
\yPx+Pxy) = -0.20345 

(PxPy) = 



TABLE IV. Correlations of the electromagnetic field vari- 
ables calculated in the fundamental state of the atom-field 
system. 



(Q+Q+) = 0.5 
(Q-Q-) = 94.05900 
(P+P+) = 0.5 
(P-P-) = 0.002657 
(Q+Q-) = 9.38113 x 10~ 10 



<(Q + P + +P+Q+)}=0 
<(Q-P-+P-Q-)} = 
((Q+P-+P-Q+)> = 
<(Q_P + +P+Q-)>=0 
(P+P-) = 4.12131 x 10~ 12 



TABLE I. Coefficients characterizing the fundamental 
state of the atom-field system calculated up to the second 
order of perturbation theory. 



an = 0.51160 012 

033 = 1 034 

ai 3 = 7.50 x 10~ 7 014 

a>23 = -5.33 x 10 -7 a 24 



= 0.78164 022 = 0.06270 

= 1.49 x 10~ 12 a 44 = 0.50751 
= 4.50 x 10~ 6 
= -7.68 x 10~ 7 
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FIG. 1. Two branches of the frequency of rotation Q(ro). 
The upper curve corresponds to the Trojan states and the 
lower curve corresponds to the anti- Trojan states; ro is mea- 
sured in units of 3600ao (ao is the atom Bohr radius). 
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FIG. 3. The boundary between the stable and unstable re- 
gions of classical equilibrium; ro is measured in units of 3600ao 
(ao is the atom Bohr radius). 
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FIG. 2. A typical classical trajectory in the rotating frame 
near the anti- Trojan equilibrium position. The motion ex- 
tends almost uniformly over the whole circle but the electron 
spends a little more time in the right half of the circle. 
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FIG. 4. Classical electron trajectory projected into the x-y 
plane; x and y are measured in units of ro. The trajectory 
started at time t = from the equilibrium position Qlq) with 
the initial momenta p = muro(0.02, k + 0.07,0.02) and is 
plotted for the time interval (1400 T, 1500 T). 
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FIG. 5. Classical motion of the electron projected into the 
phase space z-p z ; z is measured in units of ro and p z is mea- 
sured in units of muiro. The trajectory started at time t = 
from the equilibrium position ( |lq ) with the initial momenta 
p = mwro(0.02, k + 0.07,0.02) and is plotted for the time 
interval (1400 T, 1500 T). 
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FIG. 6. The energy E(n) plotted in units of hu. 




FIG. 7. The Q+-P+ cross-section of the Wigner function. 
In this (counter-rotating) mode the fluctuations of the elec- 
tromagnetic radiation are not squeezed. 




FIG. 8. The Q_-P_ cross-section of the Wigner function. 
In this (co-rotating mode) the fluctuations of the electromag- 
netic radiation are strongly squeezed. Note the change of the 
scale, as compared to Fig. Q 
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